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Notation 


R rotation matrix from a frame fixed to the 

vehicle and a global inertial reference frame. Or, if w is a vector expressed in the 
vehicle’s coordinate frame, then Rw is the same vector expressed in the global frame. 

i G , j G , k c global x-axis (North), global y-axis (East), 

global z-axis (Down) 

i A ,j A ,k A vehicle x-axis (forward), vehicle y-axis 

(right), vehicle z-axis (down) 

J inertia matrix of RLV vehicle, J e /? 3x3 , 

J = J T >0 

m vehicle mass 

g scalar, acceleration due to gravity 

q Euler angles of RLV orientation, q e R 3 (roll, 

pitch, and yaw), representing the orientation of the vehicle relative to a global inertial 
reference frame, q = (q roll , q pitch , q yaw f . 

co angular velocity of RLV expressed in 

vehicle coordinates, co<e R' 

a e R 3 , acceleration of vehicle mass center, 

expressed in vehicle coordinates, i.e., a = R T p . 

a n projection of a onto the global negative z 

unit axis -k A , i.e., the component of a , also called the “normal acceleration” of the 
vehicle 

v m vehicle Mach number (at center of mass) 

S S <e R s is the vector of control surface 

positions, for left and right flaps, left and right rudders, two inboard elevons, and two 
outboard elevons. 


^min > 


F a 

Sr 


^max e R* , s mm <0, S mwi eR\ J max > 0 , denote 

the lower and upper limits of the control surface positions S . 

j max e R" denotes an upper limit on |<j , about 


30 degrees per second for each control surface, except for the flaps with a limit of 10° V 1 . 

r a e R 3 , net aerodynamic torque exerted on 

the vehicle, expressed in the vehicle frame. 

F a & R 3 , net aerodynamic force exerted on 


the vehicle, expressed in the vehicle frame. 

drJdS , called the “torque sensitivity 

matrix.” 


A(q ) A<= R 3x3 . known partial derivative of co with 

respect q. A is nonsingular along any guidance commanded trajectory. A is a (known) 
function of q . CO = A(q)q . 
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M JA 

subscript g denotes a guidance-commanded quantity. 

For example, q roll)g is the guidance-commanded value of q rM . 

sk(v ) the 3 by 3 matrix that when post-multiplied 

by a vector w gives the cross product of v and w . Or, sk(v)w = vxw . 

p the vehicle’s linear velocity in global 

coordinates. 

Q dynamic pressure, i-/?||p||^, where p is the 

air density. 

a angle of attack 

p sideslip angle 


The Problem 


During the terminal energy management phase of flight (last of three phases) for a reusable 
launch vehicle, it is common for the controller to receive guidance commands specifying desired 
values for (i) the roll angle q roll , (ii) the acceleration a n in the body negative z direction, -k A , and 
(iii) co 3 , the projection of in onto the body- fixed axis k A , is always indicated by guidance to be 
zero. The objective of the controller is to regulate the actual values of these three quantities, i.e 
make them close to the commanded values, while maintaining system stability. 


The equations of motion are given by: 
JAq = -JAq - sk(Aq)JAq + T a 
Mq = f + T a 
q = F +M~ l T a 
ma = F a + mgR T e 3 


where 


f = -sk(Aq)JAq - JAq 
M = JA 
F = M~ l f 
e 3 = (0, 0, l) r 


( 1 ) 


The following are measured and available to the controller: S , q , in , a , a , P , v M , and Q . In 
the next section, we present a controller for this problem, including simulation results. 


Proposed Controller 


Denote the commanded orientation of the vehicle as q , whose value will be designed shortly. 
The commanded torque, t a)c , is output by the orientation controller and is chosen so as to make 
the model-based value of q satisfy: 


q + [A , ]q + [A 2 ](q~q c ) + [A 3 ] J (q - q c )dt = 0 

0 


( 2 ) 


where the [A,.] are diagonal gain matrices producing stability within the last equation. The 
values used in our simulation were: 


Aj = diag(2p,,3ju,3p), A 2 = diag(ju 2 ,3p 2 Ap 2 ), A 3 = diag(0,ju 3 ,ju 3 ) , p = 4 (3) 

giving two-percent settling times of 1 second. First, q c)l is already given as the guidance- 
command for the roll- angle. Then, the yaw-command q c)3 is chosen as q c)3 = atan2( p 2 /p t ) , which 
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would be associated with zero side-slip if p, were zero. Finally, the commanded pitch angle, 
q c)2 , is chosen as: 

9c)2 S 9 c . ( 4 ) 

where 9 C is purely model-based and approximates the value of q c)2 that would produce a n as: 

t 

a n = a n)c = a n )g ~ k i j K “ «„)* ) dt (5) 

0 

at the current mach, dynamic pressure, zero-sideslip, while q l = q c)l and q 3 = q c)3 . The value of 0 c 

is constructed as follows. Off-line, before flight, a table is constructed which inputs angle of 
attack, mach, and dynamic pressure, and outputs the corresponding a n . At the current time, t , a 
simple one-dimensional search within this table produces the value of 0 c , at virtually no 
computational cost. The second term in the expression for a n)c in (5) helps to correct for 
modeling errors. The integral gain k j was chosen as k, = 4/5 , making (5) slower than (2). 


rd 

The surface deflections S are driven by actuators with 3 order (flaps) or 
(rudders/elevons) actuator models. These actuators have either a 0.15-second 2% time-constant 
{flaps} or a 0.05-second 2% time-constant {rudders/elevons}. The commanded deflections, S c , 
are chosen as follows. Since the aerodynamic torque, r a , is not a perfectly linear function of S , 
using a delayed value of the Jacobean, Sit- A), A = 0.006.? as for one controller sample period, we 
perform a truncated {fixed number} set of damped {half- step} Newton- Rhaphson iterations, 
with the single fixed Jacobean, i.e., fixed for that time instant (or sample period). In particular, 


4 th order 



L 

r —S, ' 


<3 

II 

+ 


k 

,0.153^,0.753^ 



~ ^d)c ) J 

L 


(k=l,...,20) 


(6) 


where the subscript 1 : 8 denotes the first eight entries of the vector and where we are saturating 
at 75% of the true actuator bounds to help avoid hitting the true actuator limits. The same 
percentage is used for rate- limit avoidance. The iteration is initialized at S(t-A) . The 
simulations used 20 iterations, thus involving 20 evaluations of the look-up table for the 
aerodynamic torque (at the current v M , Q , a, and J3 = 0), 20 back substitutions, and one matrix 
factorization. The matrix G is the one associated with the stationary condition for minimizing 
the sum of the squares of the 3 i subject to achieving the specified torque: 


G = 


S(t-A) T 


(7) 


_S(t- A) 0 

The output of these Newton iterations is fed into a direct rate-limiter, with the rate-limits 

too 2 


set at 75% of v . This result is in turn fed into a filter 


s 2 +200.?+ too 2 


, the output of which is 


declared our commanded actuator deflection, S c . In all simulations, this choice of S c produced 
absolutely no saturation of the actuators’ actual bounds nor any saturation of the actuators’ 
actual rate-limits; thus, the 3 ld order (flaps) and 4 th order (rudders/elevons) dynamic actuator 
models were fully satisfied within cdl simulations, the results of which we now present in the next 
section. For possible reference, 

S min = [-15; -15; -60; -30; -30; -30; -30; -30] degrees 

3 max =(26; 26; 30; 60; 25; 30; 25; 30] degrees * ' 
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In the next section, we present numerical simulation results for a set of typical guidance 
commands. In these simulations, the guidance-commands are open-loop [which is not typical of 
the actual guidance but gives a practical test of the controller ]. 

Results 


A set of typical guidance commands (u n)d ,q roll)d ) were used to test the proposed controller. The 

nominal case’s result is shown in Figures 1-3 below. In the second case, the actual body-z air 
force is two times larger than that of the model - see Fig’s 4-6. In the third case, 1.2g {where lg 
is 9.81 m/s/s} winclgusts, each of duration 1 second, occur at times t = 75s and t = 175s — Fig’s 7- 
9. In these plots, the axes-ranges are [Os, 243s], [-70 ,20 ], [-6m/s 2 ,Sm/s 2 ] , or [-1.2 s“‘,0.4 s -1 ] , 
except for Fig’s 4, 7 and 9 with ranges [-6/n/s 2 ,10m/s 2 ] , [-10m/s 2 ,8m/s 2 ] and [-2 s“‘,4 s -1 ] . 




Figures 1-3. a n (m/s/s), roll (deg), co 3 (deg/s), dashed = actual, solid = guidance 






Figures 4-6. a n (m/s/s), roll (deg), or (deg/s), dashed = actual, solid = guidance 



Figures 7-9. a n (m/s/s), roll (deg), co 3 (deg/s), dashed = actual, solid = guidance 
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